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SUMMARY 

Stress-intensity factor equations are presented for an embedded ellipti- 
cal crack, a seml-elllptlcal surface crack, a quarter-elliptical corner crack, 
a semi-elliptical surface crack along the bore of a circular hole, and a 
quarter-elliptical corner crack at the edge of a circular hole in finite 
plates. The plates were subjected to either remote tension or bending loads. 
The stress-intensity factors used to develop these equations were obtained 
from previous three-dimensional finite-element analyses of these crack con- 
figurations. The equations give stress-intensity factors as a function of 
parametric angle, crack depth, crack length, plate thickness, and, where 
applicable, hole radius. The ratio of crack depth to plate thickness ranged 

from 0 to 1, the ratio of crack depth to crack length ranged from 0.2 to 2, 

and the ratio of hole radius to plate thickness ranged from 0.5 to 2. The 

effects of plate width on stress-intensity variations along the crack front 

were also included, but were either based on solutions of similar configura- 
tions or based on engineering estimates. 

INTRODUCTION 

In aircraft structures, fatigue failures usually occur from the initiation 
and propagation of cracks from notches or defects in the material that are 
either embedded, on the surface, or at a corner. These cracks propagate with 
elliptic or near-elliptic crack fronts. To predict crack-propagation life and 
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fracture strength, accurate stress-intensity factor solutions are needed for 
these crack configurations. But, because of the complexities of such problems, 
exact solutions are not available. Instead, investigators have had to use 
approximate analytical methods, experimental methods, or engineering estimates 
to obtain the stress-intensity factors. 

Very few exact solutions for three-dimensional cracked bodies are 
available in the literature. One of these, an elliptical crack in an infinite 
solid subjected to uniform tension, was derived by Irwin [1] using an exact 
stress analysis by Green and Sneddon [2] . Kassir and Sih [3] , Shah and 
Kobayashl [4], and Vijayakuraar and Atluri [5] have obtained closed-form 
solutions for an elliptical crack in an infinite solid subjected to non- 
uniform loadings. 

For finite bodies, all solutions have required approximate analytical 
methods. For a semi-circular surface crack in a semi-infinite solid and a 
semi-elliptical surface crack in a plate of finite thickness. Smith, Emery, 
and Kobayashl [6], and Kobayashl [7], respectively, used the alternating 
method to obtain stress-intensity factors along the crack front. Raju and 
Newman [8,9] used the finite-element method; Heliot, Labbens, and Pelllssier- 
Tanon [10] used the boundary-integral equation method; and Nishioka and 
Atluri [11] used the finite-element alternating method to obtain the same 
information. For a quarter-elliptic corner crack in a plate, Tracey [12] and 
Pickard [13] used the finite-element method; Kobayashl and Enetanya [14] used 
the alternating method. Shah [15] estimated the stress-intensity factors for 
a surface crack emanating from a circular hole. For a single corner crack 
emanating from a circular hole in a plate. Smith and Kullgren [16] used a 
finite-eleraent-alternating method to obtain the stress-intensity factors. 
Hechmer and Bloom [17] and Raju and Newman [18] used the finite-element method 
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for two symmetric corner cracks emanating from a hole in a plate. Most of 
these results were for limited ranges of parameters and were presented in the 
form of curves or tables. For ease of computation, however, results expressed 
in the form of equations are preferable. 

The present paper presents equations for the stress-intensity factors for 
a wide variety of three-dimensional crack configurations subjected to either 
uniform remote tension or bending loads as a function of parametric angle, 
crack depth, crack length, plate thickness, and hole radius (where 
applicable); for example, see Figure 1. The equations for uniform remote 
tension were obtained from Reference 19. The tension equations, however, are 
repeated here for completeness and because the correction factors for remote 
bending are modifications of the tension equations. The crack configurations 
considered, shown in Figure 2, include: an embedded elliptical crack, a semi- 

elliptical surface crack, a quarter-elliptical corner crack, a semi-elliptical 
surface crack at a circular hole, and a quarter-elliptical corner crack at a 
circular hole in finite- thickness plates. The equations were based on stress- 
intensity factors obtained from three-dimensional finite-element analyses [8, 
9, 18, and 19] that cover a wide range of configuration parameters. In some 
configurations, the range of the equation was extended by using stress- 
intensity factor solutions for a through crack in a similar configuration. In 
these equations, the ratio of crack depth to plate thickness (a/t) ranged from 
0 to 1, the ratio of crack depth to crack length (a/c) ranged from 0.2 to 2, 
and the ratio of hole radius to plate thickness (r/t) ranged from 0.5 to 2. 

The effects of plate width (b) on stress-intensity variations along the crack 
front were also Included, but were either based on solutions of similar con- 
figurations or based on engineering estimates. 
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NOMENCLATURE 


depth of crack 

width or half-width of cracked plate (see Fig. 2) 
half-length of crack 

boundary-correction factor for corner crack in a plate under tension 

boundary-correction factor for corner crack at a hole in a plate under 
tension 

boundary-correction factor for embedded crack in a plate under tension 

boundary-correction factor on stress Intensity for remote tension 

boundary-correction factor for surface crack in a plate under tension 

boundary-correction factor for surface crack at a hole in a plate 
under tension 

finite-width correction factor 

angular function derived from embedded elliptical crack solution 

curve fitting functions defined in text 

bending multiplier for corner crack in a plate 

bending multiplier for corner crack at a hole in a plate 

bending rmiltiplier on stress intensity for remote bending 

bending multiplier for surface crack in a plate 

half-length of cracked plate 

stress-intensity factor (mode I) 

applied bending moment 

curve fitting functions defined in text (i = 1, 2, or 3) 
shape factor for elliptical crack 
radius of hole 

remote bending stress on outer fiber, 3M/bt“^ 
remote uniform tension stress 

thickness or one-half plate thickness (see Fig. 2) 



function defined in text 


\ 

V Poisson's ratio ( v =* 0.3) 

^ parametric angle of ellipse, deg 


STRESS- INTENSITY EQUATIONS 

The stress-intensity factor, K, at any point along the crack front in a 
finite- thickness plate, such as that shown in Figure 1, was taken to be 

K = {S^ + HjSb)^Fj (la) 


where 



(lb) 


and 

+ (H 2 - H^) alnP ^ (Ic) 

The function Q is the shape factor for an ellipse and is given by the square 
of the complete elliptic integral of the second kind [2] . The boundary- 
correction factor, Fj, accounts for the influence of various boundaries and 
is a function of crack depth, crack length, hole radius (where applicable) , 
plate thickness, plate width, and the parametric angle of the ellipse. The 
product the corresponding bending correction. The subscript j 

denotes the crack configuration: j “ c is for a corner crack in a plate, 
j ® e is for an embedded crack in a plate, j =* s is for a surface crack in 
a plate, j *• sh is for a surface crack at a hole in a plate, and j = ch is 
for a corner crack at a hole in a plate. Functions , M 2 , M 3 , , H 2 , 

and p are defined for each appropriate configuration and loading. The 
series containing M^ is the boundary-correction factor at the maximum depth 
point. The function f^ is an angular function derived from the solution for 
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an elliptical crack In an Infinite solid. This function accounts for most of 
the angular variation in stress-intensity factors. The function f^ is a 
finite-width correction factor. Function g denotes a product of functions, 
such as gi 82 *** 8 n» used to fine-tune the equations. Functions 

and H 2 are bending multipliers obtained from bending results at (J> equal 
zero and 11 / 2 , respectively. Figure 3 shows the coordinate system used to 
define the parametric angle, 4>, for a/c less than and a/c greater than 
unity. 

Very useful empirical expressions for Q have been developed by Rawe 
(see Ref. 9). The expressions are 

1.65 

Q =» 1 + 1.464(~) for — < 1 

1.65 

Q - 1 + 1.464(-) for - > 1 

The maximum error in the stress-intensity factor caused by using these 
approximate equations for Q is about 0.13 percent for all values of a/c. 
(Rawe' 8 original equation was written in terms of a/2c). 

In the following sections, the stress-intensity factor equations for 
embedded elliptical cracks, semi-elliptical surface cracks, quarter-elliptical 
corner cracks, semi-elliptical surface cracks at a hole, and quarter- 
elliptical corner cracks at a hole in finite plates (see Fig. 2) subjected to 
either remote tension or bending loads are presented. The particular func- 
tions chosen were obtained from curve fitting to finite-element results [8, 9, 
18, and 19] by using polynomials in terms of a/c, a/t, and angular func- 
tions of (J). For cracks emanating from holes, polynomial equations in terms 
of c/r and ^ were also used. Typical results will be presented for 


(2a) 

(2b) 
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a/c » 0.2, 0.5, 1, and 2 with a/t varying from 0 to 1. Table 1 gives the 
range of applicability of <J>, a/t, a/c, r/t, and (r + c)/b for the 
proposed equations. 


Embedded Elliptical Crack 

The stress-intensity factor equation for an embedded elliptical crack In 
a finite plate. Figure 2(a), subjected to tension was obtained by fitting 
equation (1) to finite-element results In Reference 19. The results of 
Irwin [1] were used to account for the limiting behavior as a/c approaches 
zero or Infinity. The equation Is 


a „ /a a c . ^ 
^ Q ^e^c* t» b» 


(3) 


for 0 < a/c < «>, c/b < 0.5, and -it < (|) < n provided that a/t satisfies: 


- < 1.25 0.6] for 0 < - < 0 

t '^c '' c 




.2 




for 0.2 < — < « 
c 


(A) 


The function accounts for the influence of crack shape (a/c), crack size 

(a/t), finite width (c/b), and angular location ((|)), and was chosen as 


« % ^ ( 5 ) 

The terra in brackets gives the boundary-correction factors at (|> * it/2 (where 
g * ^(j) ” function f^ was taken from the exact solution for an 

embedded elliptical crack In an Infinite solid [1] and f^ is a finite-width 
correction factor. The function g Is a fine-tuning curve-fitting function. 
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For a/c < 1: 


- 1 


( 6 ) 


M 


0.05 


2 372 

0.11 + (f) 


(7) 




0.29 


0.23 + (^) 


372 


( 8 ) 


4 1/2 


( 9 ) 


and 


I J cos (|) + sin ^1 


1/4 


( 10 ) 


(Note that eq. (9) Is slightly different, and Is believed to be more accurate, 
than that given In Ref. 19.) The finite-width correction, f^^, from 
Reference 9 was 

1l/2 


( to : fa\ ^ 
w’ [®®^\,2bVt7 


( 11 ) 


for c/b < 0.5. (Note that for the embedded crack, t Is defined as one-half 
of the full plate thickness.) 


For a/c > 1: 


«1 - n/I 


( 12 ) 


8 



and 




sin <(> + cos^ 



(13) 


The functions M2, M3, g, and are the same as equations (7), (8), (9), 

and (11), respectively. 

Figure 4 shows some typical boundary-correction factors for various crack 
shapes (a/c “ 0.2, 0.5, 1, and 2) with a/t equal to 0, 0.5, 0.75, and 1. 

The correction factor, F^, is plotted against the parameter angle, (J>. At 
(1) » 0, the point on the crack front that is located at the center of the 
plate, the Influence of plate thickness is much less than at 4> » it/2, the 
point that is located closest to the plate surface. The results shown for 
a/t * 0 are the exact solutions for an elliptical crack in an infinite solid 
[1]. For a/t < 0.8, the results from the equation are within about 3 percent 
of the finite-element results. (Herein, "percent" error is defined as the 
difference between the equation and the finite-element results normalized by 
the maximum value for that particular case. This definition is necessary 
because, in some cases, the stress-intensity factor ranges from positive to 
negative along the crack front.) For a/t > 0.8, the accuracy of equation (3) 
has not been established. But its use in that range appears to be supported 
by estimates based on an embedded crack approaching a through crack (see 
Ref. 19). 

Bending equations were not developed for the embedded elliptical crack. 

Semi-Elliptical Surface Crack 

The equations for the stress-intensity factors for a semi-elliptical 
surface crack in a finite plate. Figure 2(b), subjected to remote tension and 
bending loads were obtained from Reference 9. The tension and bending 
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equations were previously fitted to finite-element results from Raju and 
Newman [8] for a/c values less than or equal to unity. Equations for 
tension and bending loads for a/c greater than unity were developed 
herein. The results of Gross and Srawley [20] for a single-edge crack 
were used to account for the limiting behavior as a/c aproaches zero. 
The equation is 


K =■ (S^ 





(U) 


for 0 < a/c < 2, c/b < 0.5, and 0 < (f> < tc, again, provided that a/t 
satisfies equation (4). The function Fg was chosen to be 


F 

8 




(15) 


For a/c < 1: 


- 1.13 - 0.09(|) 


(16) 


M 


2 


-0.54 


^ 0.89 
+ 

0.2 + - 

c 


( 17 ) 


» 0.5 


1 

0.65 + - 
c 


+ 14(1 



g - 1 + 


0 


0.1 + 0.35 



(1 


sin 



(18) 


(19) 


and is given by equation (10). The finite-width correction, f^, is 

again given by equation (11). Equations (15) through (19) were taken from 
Reference 9. (The large power in eq. (18) was needed to fit the behavior as 
a/c approaches zero.) 
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The bending multiplier, Hj, in equation (1) has the form 

+ (H^ - H^) sinP ^ (20) 

where , H 2 , and p are defined for each crack configuration considered. 

For the surface crack (j *» s), 

p » 0.2 + - + 0.6 - (21) 

c t 

- 1 - 0.34 I - 0.11 |(|) (22) 

and 

H 2 - 1 + S2(f)' (23) 

In this equation for H 2 , 

» -1.22 - 0.12 ^ (24) 

0.75 1.5 

G 22 » 0.55 - 1.05(|) +0.47(^) (25) 

Equations (21) through (25) were taken from Reference 9. 

For a/c > 1 : 

Ml “ (1 + 0.04 |) (26) 

4 

«2 “ 0.2(|) (27) 

4 

“3 “ ( 28 ) 
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(29) 


g - 1 + 


[o.l + 0.35(|)(|)^ 


(1 “ sin (J)) 


and and are given by equations (13) and (11), respectively. 


The bending multiplier for a/c >1 is also given by equation (20) where 


p * 0.2 + - + 0.6 ? 

^ at 


H = 1 + G - + G 
“l ^11 t 12 '^t'' 


(30) 

(31) 


^2 “ ^ Si t S2^t^ 


(32) 


G, , » -0.04 - 0.41 - 
11 a 


(33) 


0.75 1.5 

G^2 “ 0-55 - 1.93(^) + 


(34) 


G«, » -2.11 + 0.77 - 
21 a 


(35) 


and 


0.75 1.5 

0^2 - 0.55 - 0.72(|) + O.U(|) 


(36) 


Figures 5 and 6 show some typical boundary-correction factors for various 
surface crack shapes (a/c =■ 0.2, 0.5, 1, and 2) with a/ t equal to 0, 0.5, 
and 1 for tension and bending, respectively. For all combinations of param- 
eters investigated and a/t <0.8, equation (14) was within db5 percent of the 
finite-element results (0.2 < a/c < 2) and the single-edge crack solution 
(a/c = 0). For a/t > 0.8, the accuracy of equation (14) has not been 
established. However, its use in that range appears to be supported by 
estimates based on a surface crack approaching a through crack. 
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The use of negative stress-intensity factors in the case of bending are 
applicable only when there is sufficient tension to keep the crack surfaces 
open; that is, the total stress-intensity factor due to combined tension and 
bending must be positive. 

Quarter-Elliptical Corner Crack 

The stress-intensity factor equations for a quarter-elliptical corner 
crack in a finite plate, Figure 2(c), subjected to tension and bending loads 
were obtained by fitting equation (1) to the finite-element results in Refer- 
ence 19 for tension and the results in Table 1 for bending. The equation is 

K - (St + HcSb)f|Fc(f. (37) 


for 0.2 < a/c < 2, a/t < 1, and 0 < <{) < for ^ < 0.5, The function 

was chosen as 


F 

c 


Mi + 


M, 


(|) 


+ “3(7/ 8i 




<]> W 


(38) 


For a/c < 1: 


- 1.08 - 0.03(|) 


(39) 


M 


M 


g 


2 


3 


1 


-0.44 + 


1.06 

0.3 + - 
c 


-0.5 + 0.25(|-) + 14.8(1 - f-) 


15 


L 


1 -1 0.08 


+ 0.4(f) 


(1 - sin (j)) 


(40) 


(41) 


(42) 
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(43) 


1 + 


0.08 + 


0.15(^) 


(1 - cos <J)) 


and is given by equation (10). The finite-width correction, was 
estimated herein by using the single-edge crack tension solution given in 
Reference 21 (divided by 1.12) and was 


f « 1 - 0.2X+ 9.4X^ - 19. 4\^ + 27. IX.^ 
w 


where X * ^ V (The width correction from Ref. 21 was divided by 1.12 
because the front-face correction was already Included in eq. (38).) 
Equation (44) is restricted to c/b < 0.5. 

As a/t approaches unity, with a/c = 1 and <J) = 0, the stress- 
intensity factor equation (eq. (37)) for tension reduces to 

K * S^/tcc l.llf 
t w 

Equation (45) is within about 1 percent of the accepted solution [21] for 
c/b < 0.6. 

The bending multiplier, has the form given by equation (20). 

Functions p, , H2» ^21 given by equations (21)-(24), 

respectively, for a/c <1. The function G 22 is 

0.75 1.5 

G 22 - 0.64 - 1.05(i) +0.47(^) 


For a/c > 1: 


M 


1 



(1.08 - 0.03 |) 


M 


2 



(44) 


(45) 


(46) 


(47) 


(48) 


14 



(49) 


M 


3 



g. =* 1 + 


“ 1 + |0.08 + 


jo. 08 + 0.4(|)J 

jo. 08 + 0.15(|-)J 


(1 - sin 4>) 


(1 - cos (})) 


(50) 


(51) 


and is given by equation (13). The finite-width correction is, again, 

given by equation (44). 

The bending-correction factor is, again, given by equation (20) 

where p, , H 2 , 0^2 » and G 21 are given by equations (30)-(35), 

respectively. The function G 22 is given by 

0.75 1.5 

G.. = 0.64 - 0.72(-) + 0.14(-) (52) 

iC lU Si 


Figures 7 and 8 show some typical boundary-correction factors for corner 
cracks in plates for various crack shapes (a/c ” 0.2, 0.5, 1, and 2) with 
a/t varying from 0 to 1 for tension and bending, respectively. At a/t = 0, 
the results for tension and bending are Identical. As expected, for tension 
the effects of a/t are much larger at lower values of a/c. Again, the use 
of negative stress-intensity factors in this case of bending are applicable 
only when there is sufficient tension to keep the crack surfaces open (stress- 
intensity factor due to combined tension and bending must be positive). 

Semi-Elliptical Surface Crack at Hole 
Two symmetric surface cracks .- The stress-intensity factor equation for 
two symmetric serai-elliptical surface cracks located along the bore of a hole 
in a finite plate. Figure 2(d), subjected to tension was obtained by fitting 
equation (1) to finite-element results [19] . The equation is 
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(53) 


^ ^shC* t» t’ b» b* 


for 0.2 < a/c < 2, a/t < 1, 0.5 < r/t <2, (r + c)/b < 0.5, and 

-iu/2 < (() < n/2. (Note that here t is defined as one-half of the full plate 

thickness.) The function was chosen as 

2 


sh 


Ml + M2 (-) + M, (. 




®2 % ^ 4 » 


(54) 


For a/c < 1: 


M| “ 1 


M, 


0.05 


0.11 + (-) 


372 


(55) 

(56) 


M, 


0.29 


373 


0.23 + (^) 


(57) 




4 1/2 

(f) (2-6-2 f) 

1 + 4(^) 

''C * 


cos ((> 


(58) 


1 + 0.358\ + 1.425X,^ - 1.578X^ + 2.156X^ 
1 + 0.08\^ 


(59) 


1 + — cos(0.9 (|)) 


1 + 0.1(1 - cos (|))^(l - 


10 


(60) 


(61) 


(Note that eq. (58) is slightly different, and is believed to be more accurate, 
than that given in Ref. 19.) The function f^ is given by equation (10). 
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The finlte-width correction, f^, was taken as 


w 


sec 




2 b 


jsec 


■n:( 2 r + nc) 
4(b - c) 4- 2nc 


i-n>/2 

V t Jj 


(62) 


where n =* 1 is for a single crack and n = 2 is for two- symmetric cracks. 
This equation was chosen to account for the effects of width on stress concen- 
tration at the hole [22] and for crack eccentricity [21]. For a/c > 1: 


M 


1 



(63) 


The functions M 2 , M 3 , , g 2 , g 3 , and \ are given by equations (56) 

through (61), and the functions and f^ are given by equations (13) and 

(62), respectively. 

Estimates for a single-surface crack .- The stress-intensity factors for a 
single-surface crack located along the bore of a hole were estimated from the 
present results for two symmetric surface cracks by using a conversion factor 
developed by Shah [15]. The relationship between one- and two-surface cracks 
was given by 


(K) 


one 

crack 



(K) 


two 

cracks 


(64) 


where K for two cracks must be evaluated for an infinite plate (f^ = 1) and 
then the finite-width correction for one crack must be applied. Shah had 
assumed that the conversion factor was constant for all locations along the 
crack front, that is. Independent of the parametric angle. 

Figure 9 shows some typical boundary-correction factors for a single 
surface crack at a hole for various crack shapes (a/c “ 0.2, 0.5, 1, and 2) 
with a/t varying from 0 to 1. These results were in good agreement with 
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boundary-correction factors estimated by Shah [15] . The agreement was 
generally within about 2 percent except where the crack intersects the free 
surface (2(j>/inc = 1). Here the equation gave results that were 2 to 5 percent 
higher than those estimated by Shah. 

Stress-intensity factor equations for bending were not developed for a 
surface crack located at the center of a hole. 

Quarter-Elliptical Corner Crack at a Hole 
Two symmetric corner cracks. - The stress-intensity factor equations for 
two symmetric quarter-elliptical corner cracks at a hole in a finite plate. 
Figure 2(e), subjected to remote tension and bending loads were obtained by 
fitting to finite-element results in Reference 18. The equation is 


K 


(S, + 


«ch®b 


a r. 

F 


£. L. 

t* t* 




b» b 


(65) 


for 0.2 < a/c < 2, a/t < 1, 0.5 < r/t <2, (r + c)/b < 0.5, and 

0 < (|) < n/2. The function was chosen as 


ch 


2 4 

Mi + M2(f) +M3(f) 


®1 ®2 ®3 ^ 


w 


For a/c < 1: 


= 1.13 - 0.09 (^) 


M 


2 


“0.54 


^ 0.89 

-r 

0.2 + - 

c 


i + 14 h 

0.65 + - 
c 



( 66 ) 


(67) 

( 68 ) 


M 3 = 0.5 


(69) 



« 1 + 0.1 + 0.35(|) (1 - sin 4 >) 


1 + 0.358X+ 1.425X^ - 1.578X^ + 2.156X^ 
1 + 0 . 13 \^ 


1 + - cos( [10 


= 0.85 for tension and == 0.85 ■“ 0.25(a/t)^^^ for bending. The 
functions $3 and are given by 

, r ylf 

» (1 + 0.04 f) 1 + 0 . 1(1 - cos <t>) 0.85 + 0.15(~) 


g^ . 1 - 0 . 7(1 0 . 2)(1 -^) ( 

Functions f^ and f^ are given by equations (10) and (62), respectively. 

The bending multiplier, is given by equation (20) for a/c < 1. 

The terras p, , and H 2 are given by 


p - 0.1 + 1.3 ~ + 1.1 - - 0.7 

^ t c c '^t ^ 


2 3 

a ^ ^ r&\ 




where 


«2 ” ^^Sir-"S2(F)'^S3fr)' 


» -0.43 - 0.74 0.84 


G.o ” 1-25 - 1.19 -+ 4.39(-) 
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(80) 


2 

Gi 3 “ -1.94 + 4.22 5.5l(^) 


2 

G 21 = -1.5 - 0.04 - 1*73(^) 


2 

G 22 =* 1.71 - 3.17 6.84(^) 


2 

G 23 =» -1.28 + 2.71 5.22(^) 


For a/c > 1 : 


Ml - /I (1 + 0.04 f) 
4 

M 2 = 0.2 (|) 

4 

M3 » -0.11(f) 


81 


1 + 


0.1 + 0.35(^)(|-) 


(1 - sin (|>) 


(81) 

(82) 

(83) 

(84) 

(85) 

( 86 ) 

(87) 


Functions g 2 and X are given by equations (71) and (72). Function g 3 is 
given by 


83 . (1.13 - 0.09 |) 


1 + 0 . 1(1 ~ cos (j))' 


0.85 + 


0- 15(f) 


1/4 


( 88 ) 


and g^ *■ 1. The functions and are, again, given by equations (13) 

and (62), respectively. 

Again, the bending-correction factor, is given by equation (20). 

The function p is given by equation (30) for a/c > 1. The H-functions are 
given by equations (75) and (76) where 


20 



(89) 


“ -2.07 + 0.06 j 


Gt 2 “ 4,35 + 0.16 - (90) 
G ^3 - -2.93 - 0.3 ~ (91) 
G 21 “ “3.64 + 0.37 ~ (92) 
G 22 “ 5.87 - 0.49 ~ (93) 
G 23 = -4.32 + 0.53 ~ (94) 


Eatlmates for a single-corner crack .- The stress-intensity factors for a 
single-corner crack at a hole were estimated from the present results for two- 
symmetric corner cracks by using the Shah-conversion factor (eq. (64)). Raju 
and Newman [18] have evaluated the use of the conversion factor for some 
corner-crack-at-a-hole configurations. The stress-intensity factor obtained 
using the conversion factor were in good agreement with the results from Smith 
and Kullgren [16] for a single-corner crack at a hole. 

Figures 10 and 11 show some typical boundary-correction factors for a 
single corner crack at a hole for various a/c and a/t ratios for tension 
and bending, respectively. Again, the use of negative stress-intensity 
factors in the case of bending are applicable only when there is sufficient 
tension to make the total stress-intensity factor, due to combined tension and 
bending, positive. 
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CONCLUDING REMARKS 


Stress-intensity factors from three-dimensional finite-element analyses 
were used to develop stress-intensity factor equations for a wide variety of 
crack configurations subjected to either remote uniform tension or bending 
loads. The following configurations were included: an embedded elliptical 

crack, a serai-elliptical surface crack, a quarter-elliptical corner crack, a 
semi-elliptical surface crack along the bore of a hole, and a quarter- 
elliptical corner crack at the edge of a hole in finite plates. The equations 
cover a wide range of configuration parameters. The ratio of crack depth to 
plate thickness (a/t) ranged from 0 to 1, the ratio of crack depth to crack 
length (a/c) ranged from 0.2 to 2, and the ratio of hole radius to plate 
thickness (r/t) ranged from 0.5 to 2 (where applicable). The effects of plate 
width (b) on stress-intensity variations along the crack front were also 
included, but were based on engineering estimates. 

For all configurations for which ratios of crack depth to plate thickness 
do not exceed 0.8, the equations are generally within 5 percent of the finite- 
element results, except where the crack front intersects a free surface. Here 
the proposed equations give higher stress-intensity factors than the finite- 
element results, but these higher values probably represent the limiting 
behavior as the mesh is refined near the free surface. For ratios greater 
than 0.8, no solutions are available for direct comparison; however, the 
equations appear reasonable on the basis of engineering estimates. 

The stress-intensity factor equations presented herein should be useful 
for correlating and predicting fatigue-crack-growth rates as well as in 
computing fracture toughness and fracture loads for these types of crack 
configurations. 
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Table 1. Range of applicability for stress-intensity factor equations. 


a^ 


Configuration 

Equation 



a/ 1 

a/c 

r/t 

(r + c)/b 

Embedded Crack in Plate^®^ 

(3) 

-7T 

to ir 

(b) 

0 to » 

• • • 

<0,5^^^ 

Surface Crack in Plate 

(14) 

0 

to IT 

(b) 

0 to 2 

« ♦ • 

<0.5^^^ 

Corner Crack in Plate 

(37) 

0 


<1 

0.2 to 2 

... 

<0.5^^^ 

Surface Crack at Hole^^*^^ 

(53) 

TT 

“ 2 

tol 

<1 

0.2 to 2 

0.5 to 2 

<0.5 

Corner Crack at Hole^^^ 

(65) 

0 

tof 

(e) 

0.2 to 2 

0.5 to 2 

<0.5 


(a) 

(b) 

(c) 

(d) 

(e) 


Equations for bending were not developed for this case, 
a/t < 1.25 (a/c + 0.6) for 0 < a/c < 0.2 and a/t < 1 for 
r = 0. 

One or two-s 5 nnnietric cracks. 

a/t < 1 for remote tension and a/t < 0.8 for remote bending 


a/ 


c 


> 0 . 2 . 



Table 2. Boundary-correction factors, F^H^, for quarter-elliptical 


corner 
= 0 

crack in 

.3; F H = 
c c 

a plate subjected 
K/(S^vW^ )). 

to bending 





a/t 


a/c 

2(j)/7r 

0.2 

0.5 

0.8 


0 

0.522 

0.609 

0.779 


0.25 

0.669 

0.702 

0.808 

0.2 

0.5 

0.801 

0.746 

0.716 


0.75 

0.868 

0.746 

0.577 


1.0 

0.876 

0.750 

0.604 


0 

0.740 

0.799 

0.904 


0.25 

0.724 

0.690 

0.670 

0.4 

0.5 

0.785 

0.632 

0.451 


0.75 

0.826 

0.583 

0.272 


1.0 

0.846 

0.569 

0.262 


0 

1.084 

1.046 

1.027 


0.25 

0.934 

0.770 

0.604 

1.0 

0.5 

0.838 

0.547 

0.237 


0.75 

0.798 

0.417 

0.011 


1.0 

0.839 

0.407 

-0.032 


0 

0.932 

0.811 

0.734 


0.25 

0.851 

0.623 

0.442 

2.0 

0.5 

0.761 

0.413 

0.105 


0.75 

0.700 

0.268 

-0.131 


1.0 

0.677 

0.215 

-0.206 
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Fig. 1 — Corner cracks at the edge of a hole in a 
finite plate subjected to remote tension 

and bending. 
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- ^ 

(c) Corner crack 


H — 2r-H c 

(e) Corner crack at hole 


Fig. 2— Embedded-, surface-, and corner-crack configurations (all cracks have 
elliptical fronts). 




















20/ 7T 20/ 7T 


(c) a/C = 1 . (d) a/C = 2 . 


Fig. 7 — Typical boundary-correction factors for a corner crack in a 
plate subjected to remote tension (c/b = 0). 
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0 


.5 

2(p/n 

(c) a/c = 1. 


1 


0 


.5 

2<p/ 

(d) a/c = 


Typical boundary-correction factors for a single si 
at the center of a circular hole In a plate subject 
tension (r/t = 1; r/b = 0), 



u/t 


U./J. 


(b) a/c = 0.5. 



Fig. 10— Typical boundary-correction factors for a single corner crack 

at the edge of a circular hole In a plate subjected to remote 
tension (r/t - 1: r/b = OK 
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